Multi-Particle Collision Dynamics -- a Particle-Based Mesoscale Simulation Approach to the Hydrodynamics of Complex Fluids

In this review, we describe and analyze a mesoscale simulation method for fluid flow, which was introduced by Malevanets and Kapral in 1999, and is now called multi-particle collision dynamics (MPC) or stochastic rotation dynamics (SRD). The method consists of alternating streaming and collision steps in an ensemble of point particles. The multi-particle collisions are performed by grouping particles in collision cells, and mass, momentum, and energy are locally conserved. This simulation technique captures both full hydrodynamic interactions and thermal fluctuations. The first part of the review begins with a description of several widely used MPC algorithms and then discusses important features of the original SRD algorithm and frequently used variations. Two complementary approaches for deriving the hydrodynamic equations and evaluating the transport coefficients are reviewed. It is then shown how MPC algorithms can be generalized to model non-ideal fluids, and binary mixtures with a consolute point. The importance of angular-momentum conservation for systems like phase-separated liquids with different viscosities is discussed. The second part of the review describes a number of recent applications of MPC algorithms to study colloid and polymer dynamics, the behavior of vesicles and cells in hydrodynamic flows, and the dynamics of viscoelastic fluids

Lattice Boltzmann simulations of soft matter systems

This article concerns numerical simulations of the dynamics of particles immersed in a continuum solvent. As prototypical systems, we consider colloidal dispersions of spherical particles and solutions of uncharged polymers. After a brief explanation of the concept of hydrodynamic interactions, we give a general overview over the various simulation methods that have been developed to cope with the resulting computational problems. We then focus on the approach we have developed, which couples a system of particles to a lattice Boltzmann model representing the solvent degrees of freedom. The standard D3Q19 lattice Boltzmann model is derived and explained in depth, followed by a detailed discussion of complementary methods for the coupling of solvent and solute. Colloidal dispersions are best described in terms of extended particles with appropriate boundary conditions at the surfaces, while particles with internal degrees of freedom are easier to simulate as an arrangement of mass points with frictional coupling to the solvent. In both cases, particular care has been taken to simulate thermal fluctuations in a consistent way. The usefulness of this methodology is illustrated by studies from our own research, where the dynamics of colloidal and polymeric systems has been investigated in both equilibrium and nonequilibrium situations.Comment: Review article, submitted to Advances in Polymer Science. 16 figures, 76 page

Shear-induced migration in flowing polymer solutions: Simulation of long-chain deoxyribose nucleic acid in microchannels

This material is presented to ensure timely dissemination of scholarly and technical work. Copyright and all rights therein are retained by authors or by other copyright holders. All persons copying this information are expected to adhere to the terms and constraints invoked by each author's copyright. In most cases, these works may not be reposted without the explicit permission of the copyright holder

Recirculating Flows Involving Short Fiber Suspensions: Numerical Difficulties and Efficient Advanced Micro-Macro Solvers

Numerical modelling of non-Newtonian flows usually involves the coupling between equations of motion characterized by an elliptic character, and the fluid constitutive equation, which defines an advection problem linked to the fluid history. There are different numerical techniques to treat the hyperbolic advection equations. In non-recirculating flows, Eulerian discretizations can give a convergent solution within a short computing time. However, the existence of steady recirculating flow areas induces additional difficulties. Actually, in these flows neither boundary conditions nor initial conditions are known. In this paper we compares different advanced strategies (some of them recently proposed and extended here for addressing complex flows) when they are applied to the solution of the kinetic theory description of a short fiber suspension fluid flows

Mobility gradient induces cross-streamline migration of semiflexible polymers

Many aspects of modern material science and biology rely on the strategic manipulation and understanding of polymer dynamics in confining micro- and nanoflow. We directly observe and analyze nonequilibrium structural and dynamic properties of individual semiflexible actin filaments in pressure-driven microfluidic channel flow using fluorescence microscopy. Different conformational shapes, such as filaments fluctuating in an elongated manner, parabolically bent, as well as tumbling, are identified. With increasing flow velocity, a strong center-of-mass migration toward the channel walls is observed. This significant migration effect can be explained by a shear rate dependent spatial diffusivity due to a gradient in chain mobility of the semiflexible polymers